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Abstract 

We present the first numerical application of a method that we have recently proposed to 
solve the Non Perturbative Renormalization Group equations and obtain the n-point functions 
for arbitrary external momenta. This method leads to flow equations for the n-point functions 
which are also differential equations with respect to a constant background field. This makes 
them, a priori, difficult to solve. However, we demonstrate in this paper that, within a simple 
approximation which turns out to be quite accurate, the solution of these flow equations is not 
more complicated than that of the flow equations obtained in the derivative expansion. Thus, 
with a numerical effort comparable to that involved in the derivative expansion, we can get 
the full momentum dependence of the n-point functions. The method is applied, in its leading 
order, to the calculation of the self-energy in a 3-dimensional scalar field theory, at criticality. 
Accurate results are obtained over the entire range of momenta. 
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I. INTRODUCTION 

The non perturbative renormalization group (NPRG) |2|, |3j, |4|, |5( stands out as a 
very promising formalism to address non perturbative problems, i.e., problems in which 
the absence of a small parameter prevents one to build a solution in terms of a system- 
atic expansion. It leads to exact flow equations which are difficult to solve in general, 
but which offer the possibility for new approximation schemes. When only correlation 
functions at small momenta are needed, as is the case for instance in the calculation of 
critical exponents, a general approximation method to solve the infinite hierarchy of the 
NPRG equations has been developped 0, Q| • This method, which can be systemati- 
cally improved, is based on a derivative expansion of the the effective action. It has been 
applied successfully to a variety of physical problems, in condensed matter, particle or 
nuclear physics (for reviews, see e.g. p, Q]). However, in many situations, this is not 
enough: in order to calculate the quantities of physical interest, the knowledge of the full 
momentum dependence of the correlation functions is mandatory. Many efforts to get this 



information from the flow equations, involve truncations of various kinds 
early suggestion by Weinberg (see however till ) • 
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following an 
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'he present paper explores the ap plic ability of the strategy that we proposed recently 
in |l4 1 , following our previous works jisi in which we presented a scheme to obtain 
the momentum dep endence of n-point functions from the flow equations. The strategy 



He 
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put forward in [1^ is based on the fact that the internal momentum q in the integrals 
that determine the flow of the n-point functions is bounded by the regulator introduced 
by the NPRG. Since this regulator also guarantees that the vertex functions are smooth 
functions of the momenta, these can be expanded in powers of q 2 /^ 2 , k being the cut-off 
scale in the regulator. The "leading order" (LO) of the approximation scheme proposed 
in [l^l simply consists in keeping the lowest order of this expansion, i.e., in setting q = 
in the vertices. Doing so, and working in a constant external field, it is possible to relate 
to each other the various n-point functions that appear in a given flow equation through 
derivatives with respect to the external field, thereby closing the hierarchy of equations. 

In we showed that the method reproduces perturbative results, at any desired 
order. Furthermore, we also showed that the LO is exact in the large N limit of the O(N) 
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scalar model. Finally, one expects the method to provide, at each order of the expansion, 
results as good as those of the derivative expansion in the domain where the derivative 
expansion is valid. 

The price to pay is that the resulting equations are also differential equations with 
respect to a uniform background field, with integral kernels that involve the solution itself. 
These integro-differential equations are a priori difficult to solve. The aim of this paper 
is to demonstrate that they can indeed be solved, with a numerical effort comparable to 
that involved in solving the flow equations that result from the derivative expansion, and 
to present a first application to the study of the 2-point correlation function of the scalar 
model, in the LO of the approximation scheme. 

The outline of the paper is as follows. In section |H] we briefly recall some basic 
features of the NPRG and the essence of our approximation scheme in the case of a scalar 
field theory. In section ITTT1 we analyze the structure of the flow equation for the 2-point 
correlation function and describe the strategy that we used to solve it. In section ITVl we 
present numerical results for the self-energy of the scalar field, at criticality and in d = 3. 
The appendices gather technical material. 

II. THE METHOD 

We consider a scalar field theory with the classical action 

S = jd d x {i(^(x)) 2 + ^ 2 (x) + ^ 4 (x)J . (1) 

The NPRG constructs a family of effective actions, (with <fr the expectation value of 

the field in the presence of external sources), in which the magnitude of long wavelength 
fluctuations are controlled by an infrared regulator depending on a continuous parameter 
k. The effective action interpolates between the classical action obtained for k = A 
(with A the microscopic scale at which fluctuations are essentially suppressed), and the 
full effective action obtained when k — > 0, i.e., when all fluctuations are taken into account 
< S ee e., Q). It is understood that the va.ne 8 of the palters r and u of the classical 
action as well as the field normalisation, are fixed at the microscopic scale A. One 
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can write for T K [(j>] an exact flow equation 

1 f d d q 



(27T) 



(9 KJ R K (g))[ri 2 ) + J R K ] 



-i 
q-q 



(2) 



where r« is the second derivative of T K with respect to 0, and R K denotes a family 



of "cut-off functions" depending on k. 
abundantly discussed in the literature jl7 



rere is a large freedom in the choice of R K (q), 



18, 
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To be specific, in the present paper, 



we shall use for R K (q) the following function 

R K {q) = Z.{K 2 -q 2 )e{K 2 -q 2 ) 



(3) 



where Z K is a function of k specified in the next section (see eq. (JUJ)). 

By deriving eq. (J2J) with respect to 0, and then letting the field be constant, one gets 
the flow equation for the n-point functions r^ n ^ in a constant background field 0. More 
precisely, taking into account momentum conservation, one defines: 



(2n) d tfM ( Pl + ■ • • + p n ) rW( Pl ,...,p n ; 



S n T K 



54>{xi) . ..5(j){x n ) 



(4) 



Then, the equation for the 2-point function reads: 

d d q 



P- 



(2vr) ( 



{d K R K {q)) {G K (q; <P)^\p, q, ~P ~ q; 0) 

xG K (q + p; <P)Tf{-p,p + q, -q; 0)G K (g; 0) 



(5) 



where 



G-\ q]< P)^T^\q,-q-A) + RM, (6) 
and in eq. (JSJ) we have used the simplified notation Y K (q; 0) for T K (q, — q; 0), a notation 
that will be used throughout. 

In general, the flow equation for a given n-point function involves the m-point func- 
tions with m = n + 1 and m = n + 2. Thus, the flow equations for the n-point functions 
do not close, but constitute an infinite hierarchy of coupled equations; this makes them 
difficult to solve. 

In 3] we proposed a method to solve this infinite hierarchy. It exploits the smooth- 
ness of the regularized n-point functions at small momenta, and the fact that the loop 



momentum q in the right hand side of the flow equations (such as eq. (J2J) or eq. (|5|)) is 
limited to q ^ k by the presence of the regulator R R (q). The leading order (LO) of the 
method presented in thus consists in setting q = in the n-point functions in the 
r.h.s. of the flow equations, for instance 

rt n) (Pl,P2, ...,p n -l + q,Pn~q) ~ ll n) (pi,p2, p n -l, Pr>) ■ (7) 

Once this approximation is made, some momenta in some of the n-point functions vanish, 
and the corresponding n-point functions can be obtained as the derivatives of m-point 
functions (m < n) with respect to a constant background field, thereby allowing us to 
close the hierarchy of equations. 

Specifically, in eq. (0) for the 2-point function, the 3- and 4-point functions in the 
r.h.s. will contain respectively one and two vanishing momenta after we set q = 0. These 
can be related to the following derivatives of the 2-point function: 

One then arrives at a closed equation for r» (p; p) (with p = 2 /2): 

p r(2)^ \ 7 (% n /dil 2) (p;p)\ l r (2)/ ,d 2 T^\p;p) 
Kd K r^{p;p) = jy{p;K;p)[ — I --^ ; (/?;p) — , (9) 

where 

4 n \p;^p)^ f ^A d M ( l))G K (p + q; P )G^- 1 \q-,p), (10) 

and 

#W) = J -0^K(d K R K (q))G n K (q; p). (11) 

Note that j£°(p = 0; «; p) = ^(k; p). 

At this point we note that the n-point functions at zero external momenta can all be 
considered as derivatives of a single function, the effective potential V K (p). For instance, 

rf(P = 0;p) = ^f (12) 
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The effective potential satisfies a flow equation which can be deduced from that for the 
effective action, eq. (J2J), when restricted to constant fields. It reads 

K8«V K {p) = \j j^K(d K R K (q))G K (q;p). (13) 

The second derivative of this equation with respect to the background field yields a flow 
equation for T K (p = 0; p). Now, this equation does not coincide with eq. © in which we 
set p = 0: indeed, in contrast to eq. (JSJ), the vertices in the equation deduced from eq. (fT5Jl 
keep their g-dependence (q being the loop momentum in eq. (|13)0. There is therefore an 
apparent inconsistency in our approximation scheme, that is however easily resolved by 
treating separately the zero momentum (p = 0) and the non-zero momentum (p ^ 0) 
sectors. In fact, in doing so, we get more accuracy in the sector p = than in the sector 
p^O. 

Let us then write : 

ri 2) (p;p)=p 2 + ^ + S K (p;p), (14) 

where 

£ K (p; p) = rW(p; p)-p 2 - Vf{p = 0; p). (15) 
We shall refer to £ K (p; p) as the self-energy (although it differs from the usual self-energy 

(2) 

by the subtraction of the momentum independent contribution T K (p = 0;p)). By defi- 

(2) 

nition, £ K (p = 0; p) = 0, and at criticality, T K: L (p = 0; p) = 0. We then proceed with 
separate approximations in the two sectors with p = and p ^ 0. 

In the sector p ^ 0, it is £ K (p; p), rather than T K (p; p) which satisfies the approximate 
eq. (jHJ) (strictly speaking, eq. (jHJ) to which one subtracts the same equation with p = 0): 

-[p-0]. (16) 

Eq. (jlfjjl is the flow equation for the momentum dependent part of the 2-point function 
at LO of our approximation scheme. It is a partial differential equation with respect to 
the two real variables, k and p, with the momentum p playing the role of a parameter. It 
is to be integrated from k — A, with initial condition £a(p; P) — ( see ec L s - © an d (US))), 
to k = where it yields the physical self-energy E(p; p) = S K=0 (p; p). 
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ndjl(p]p) 



jf\p- K-p) 



arffep) 



2 Id 



(2), 



K-p) 



a 2 rl 2) (p;p) 



Eq. (JBj) is to be solved together with the equation in the sector p — 0, i.e., with 
the equation for the effective potential, with initial condition V\(p) = rp + (w/6)p 2 (see 
eq. (JU)). In eq. (fT3j) for V K (p), we use the propagator (fTlj) in which £ K (p;p) is solution 
of eq. (|TH|) and I« (g = 0; p) is determined self-consistently from the effective potential, 
using eq. fH2|) . 

It is not difficult to verify that (in the perturbative regime) this scheme has 2-loop ac- 
curacy for the effective potential, and only one-loop accuracy for the self-energy. Besides, 
in the low momentum region it is as accurate as the derivative expansion at next-to-leading 
order. 

III. ANALYSIS OF THE FLOW EQUATION 

There are two features of eq. |T6|) that make it a priori difficult to solve. First, 
the two functions J d (p; k; p) and I d {k\ p), are functional of the solution T K (p; p) (see 
eq. (P)). Second the different values of p are coupled through the propagator G K (p + q) 
entering the calculation of J^\p;K;p). In principle, one should therefore solve eq. (fTH|) 
self-consistently, and simultaneously for all values of p. However, in this section, we shall 
show that it is possible to make an accurate calculation of jf\p] k; p) and I d {K\ p) using 
approximate propagators. This yields an approximate version of eq. (|T6|) that can be 
solved for each given value of p. The validity of this approximation will be checked in the 
next section. 

Consider first the function I^\k; p), which does not depend on p. The smoothness of 
the n-point functions and the fact that q < k, suggest to perform in the propagators of the 
right-hand-side of eq. (fTTj) an approximation similar to that applied to the other n-point 
functions, i.e., set q = 0. However, in order to maintain the exact one-loop properties of 
the flow equations, one cannot simply set q = in the propagators: rather, one needs to 
keep a momentum dependence close to that of the free propagators. Thus, we shall use 
for the propagators entering the calculation of I d n \n] p) the following approximate form 

G-\q; p) « Z K q 2 + V®{q = 0; p) + R K (q), (17) 
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where 



dq 2 



q=o,p=po 



As well known ja], and will be verified in app. A, dT^(q; p)/dq 2 \^_ Q depends weakly on 
p. Accordingly, one expects Z K to depend weakly on the value chosen for p . As will be 
seen in app. A, the choice po = is here the simplest. With the propagator ()17j) . and the 
function Q for R K (q) one can calculate I d n \n; p) analytically: 

In this expression, 

VK = -ndJnZ K (20) 

is the running anomalous dimension and 

»2/ ,_ ri 2) (g = 0;p) 

"W) = , (21) 

is a dimensionless, field-dependent, effective mass. is a number resulting from angular 
integration, Kj 1 = d 2 d ~ l ix d / 2 T(d/2) (e.g., K 3 = 1/(6tt 2 )). Notice that, for d > 2, 
I d [K] p) — > when k — >■ 0. 

We shall calculate the function J d (p, K] p) in a similar way, arguing that in this 
calculation one can assume p ^ k : the propagator G K (p + q; p) in eq. (fTUj) is small as 
soon as p/n is large, and one can indeed verify that the function J d (p; k; p) vanishes 
approximately as K 2 /p 2 for large values of p/n (see the explicit expressions f|B3j) and ()B5j) 
given in app. B). Thus, in the region where J d (p; k; p) has a significant value, one can 
use for G K {p + q; p) an expression similar to (JT7J), namely 

^(p + q; p) « Z K (p + g) 2 + if >(0; p) + i? K (p + g). (22) 

One can then calculate the function J d 3 \p; k; p) analytically (in d — 3). The resulting 
expression is more complicated than that of I d 2 \n; p), eq. (JTHJ). It is given in app. |B] (see 
also [3|). Observe that the regulator in eq. (jHJ) is not analytic at q ~ k. This generates 
non analyticities in J d (p; k; p); but these occur only in the third derivative with respect 



9 



to p, at p = and at p = Ik (cf. the odd powers of p in eqs. flB3HB6jO . and they play no 
role at the present level of approximation. 

(2) (3) 

With the approximations just discussed, 7^ (k\ p) and J d (p; k\ p) depend only on 
quantities that enter the flow equations at p — 0, namely mf.(p) and Z K (or rj K ). As we 

discuss in app. 1X1 these quantities can be obtained from a modified version of the Local 

1 

Potential approximation |7|] that we call the LPA'. The strategy to solve eq. consists 
then in two steps: one first solves the LPA' to get rn? K {p) and r] K ; then, for each value of 
p, one solves eq. (fTBj) with the kernels if (k; p) and jf (p\ k; p) that are calculated with 
rh^(p), Z K and rj K determined from the LPA'. 

Note that, generally, the flow of £ gets strongly suppressed below some non vanishing 

(2) 

value of k. In d = 3, this can be inferred from the properties of the functions 1% (k; p) 
and Jg (p; k; p) discussed above, and it will be verified explicitly on the numerical results 
presented in the next section (see Fig. P). In fact, the flow of £ receives two contribu- 
tions: the first involves the external momentum p ^ and is suppressed when k ^ p 

(3) (2) 

(J3 (p; k\ p) vanishes rapidly when k becomes smaller than p, while Jg (k; p) ~ k); the 
other contribution is independent of p and, at criticality, is suppressed for n ^ k c ~ u/10 
(see app. A, and in particular fig. |HJ). Accordingly, one expects the flow to stop when k 
reaches the smallest of k c and p. 

The function (p~, p) exhibits a simple scaling behavior. Consider for simplicity the 
zero field case p = 0, and the ratio 



At criticality, and in the scaling regime where p, k u, we expect / to become indepen- 
dent of u, and therefore a function of p/n only. As will be shown in the next section, the 
solution of eq. (j!6|) verifies this property. Note that this scaling behavior is reproduced 
only when including a renormalization factor Z K whose flow is determined consistently 
from that of dH K /dp 2 for p < k, as obtained from eq. (J16)) . This calculation of Z K is 
explained in app. A. We have tested the consequence of setting Z K = 1 in the propagators 
(JT7J) and (jUJ), corresponding to the Local Potential approximation (as opposed to the 
LPA'). Doing so does not alter the self-energy in any significant way when p ^ u, but in 
the IR regime, the scaling behavior is only approximate. 




(23) 
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FIG. 1: Left: p = 0)/it 2 as a function of n/u for various values of p 3> u. The flow stops 

at k 5i k c ~ it/10. Right: p = 0)/E re= o(p; p = 0) as a function of n/p for various values of 

p <C it. The flow stops at k ^ p/10. 

IV. NUMERICAL RESULTS AND DISCUSSION 

We now turn to the numerical solution of the flow equation for S K (p; p), at d — 3 and 
at criticality. Our goal is to assess the quality of the approximation scheme, and there 
are two aspects that we shall examine. First, since the strategy described in the previous 
section provides only an approximate solution to eq. (jlfij) . we shall estimate by how much 
this approximate solution differs from the exact solution of this equation. Second, since 
eq. (fhoj) itself is only the LO approximation of the method described in [14j, we shall 
compare our results with known ones concerning the self-energy of the scalar model at 
criticality. 

Let us start by considering general properties of the flow, and verify in particular that 
it essentially stops for a small value of k. Fig. [T] displays the self-energy, p = 0) as 

a function of k/u, for different values of p. Calculations are made for u/A = 3.54 x 10~ 4 
(this value is small enough to guarantee that the results are independent of A). The 
left panel of fig. ^ shows the flow of p = 0) for values of p in the UV regime, i.e., 
p ^> n c ~ u/10; for all the considered values of p the flow stops at k c . The right panel of 
fig. presents the flow of the self-energy when p is in the IR regime, i.e., when p k c . In 
this case, we have divided £ K (p; p = 0) by its physical value E(p; p = 0) = S K=0 (j9; p = 0), 
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FIG. 2: T,(p)/u 2 , in d = 3, at criticality and zero external field, as a function oip/u. 

in order to make it more obvious that the flow only stops when k ^ p/10. 

We now turn to the physical self-energy = £ K=0 Q»; p = 0) in vanishing external 
field, displayed in fig. El as a function of p/u, and discuss its behavior in the various 
momentum regions: p > ti, p < ti, p ~ k c ~ ?i/10. We have checked that the curve 
in fig. El i-e., (l/u 2 )S(p/w), is "universal", i.e., independent of u and A, provided u/A is 
small enough. 

In the perturbative regime (p ^> u), one expects ~ (u 2 /96n 2 ) log(p/u). In app.lUl 
we show that the analytical solution of eq. ()16|) preserves this behavior, although the 
coefficient in front of the logarithm is (u 2 /9ir 4 ), 8% larger (the LO approximation does not 
include all the 2-loop perturbative diagrams exactly). Our approximate numerical solution 
reproduces this result. As explained in Q|, at the NLO of our approximation scheme, 
which is beyond the scope of the present paper, the contribution of the 2-loops diagrams 
would be exactly included and the correct prefactor (u 2 /96tt 2 ) would be recovered. 

In the IR region (p <ti u) we expect the self-energy to behave as 

p 2 + Z(p) = Ap 2 -v\ (24) 
where rj* is the anomalous dimension. By analyzing the small momentum behavior of 
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FIG. 3: The ratio (p 2 + T, K (p))/(Z K p 2 ) as a functio of p/k. Right: the same ratio divided by 
(p/«)-"*. 

we get numerically 77* = 0.05218. An alternative way to determine i]* is to extract 
it from the k dependence of Z K (see eq. (I20j) ). As recalled in app. A, in the critical regime, 
i.e., when n * , with r]* = 0.05220 the fixed point value of r] K (see fig. EJ). 

It is also shown in fig. E]m app. A that the quantity T^\p = 0; p = 0)/(k 2 Z k ) goes to 
a fixed point, which confirms the behavior of T^\p = 0; p = 0) ~ k 2 ~ v expected in the 
scaling regime. 

We have performed a more stringent test of scaling by studying the function 
{p 2 + S K (p))/(Z K p 2 ). This function is displayed in fig.Elas afunction oip/K. By definition 
of Z K (see eq. fITHj) ). when k is kept fixed and p — > 0, this function goes to one. Further- 
more, as explained before, in the scaling regime p,n <C u, one expects this function to 
depend on p/ k only, which is indeed well verified, as can be seen on the left panel of fig. El 
it is only for values of p which are not small enough (p/u = 10~ 2 , corresponding to the 
dashed line) that violation of this scaling starts to become significant. Moreover, as can 
be seen on the figure, p 2 + S K (p) is well approximated by Z K p 2 for all p ^ k. In the right 
panel of fig. El we have plotted the ratio (p 2 + T, K (p))/(Z K p 2 ) divided by (p/n)' 71 * . Recall 
that when k < p < k c , one expects p 2 + S K (p) ~ p 2 ~ v * , while Z K ~ k~ v * when k ^ k c . 
Therefore when 1 < p/k < n c / k, one expects (p 2 + £ K (p))/ \Z K p 2 ) ~ (p/k)~ v * , so that 
the quantity which is plotted should be constant. As seen in the right panel of fig. El this 
is indeed the case for the value rj* = 0.05219, which confirms the coherence of the whole 
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calculation. 

Our estimate for the anomalous dimension, 77* rs 0.052 is to be compared with the 
results 77* — 0, rj* — 0.044 and 



NLO and NNLO, respectively 



12. 



0.033 obtained with the derivative expansion at LO, 
, and 77* = 0.0335 ± 0.0025 from the resummed 7 
loop calculation of ref. [24|. Thus, the LO of our approximation scheme yields a result 
slightly larger that the NLO of the derivative expansion. The value of 77* obtained here is 
also slightly larger than that obtained in ^| using a different version of the LPA' than 
that used here. In fact, the deviation of the present estimate of 77* from the value 0.044 
obtained with the derivative expansion in next-to-leading order can be taken as a measure 
of the error introduced, in the scaling regime, by our use of the LPA' in our approximate 
solution of eq. (fTHjl: as already mentioned, if we had solved eq. (fTfij) exactly, one should 
have obtained essentially the same value as in the derivative expansion at NLO. 

We turn now to the intermediate momentum region, which we shall probe with a 
quantity which is very sensitive to the cross-over between the two regimes just studied: 

^W^f^-iY (25, 



As shown for instance in 



(2tt) 3 \p 2 + Z(p) p 2 / 
2l|, the integrand in eq. ()25j) is peaked at p ~ k c (in fact it 

takes significant values only in the region 10~ 3 ^ p/u ^ 10). This quantity has been much 

studied recently for a scalar model with 0(2) symmetry because it determines then the 

shift of the critical temperature of the weakly repulsive B ose g as [25| . For the simple scalar 

model studied in this paper, lattice calculations measure 12611 A((p 2 )/u = —(4.95 ±0.41) x 

10~ 4 while the "7 loop" resummed calculation of ref. [27[ yields A{<p 2 )/u = —(4.86 ± 

0.45) x 10~ 4 . With the present numerical solution, one gets A(<fr 2 )/u = —5.39 x 10~ 4 . 

This is only slightly larger than the value A((t> 2 ) /u = —5.03 x 10~ 4 obtained in the next- 



to-leading order of the scheme presented in |l5l . Il6 1 . 



We conclude that with the LO of the present approximation scheme, we obtain an 
accurate description of the self-energy in the entire range of momenta. Since we have 
solved only appproximately eq. it remains to study by how much the solution that 
we have obtained differs from the exact solution of eq. (|TH|) . We have already indications 
about the accuracy of the approximation both in the UV and in the IR. In the UV, we 
reproduce the expected result (which differs by 8% from the exact 2-loop result). We 
also loose the 2-loop accuracy with which the efffective potential could be obtained in 
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k/u=3 10" 
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(2) 

FIG. 4: Left: The ratio of the function Ig (k; p) calculated with the obtained numerical 
propagator and with the approximate LPA' propagator (as explained in the text), as a function of 
p, for different values of k/u. Right: The function 1^ (ft; p) as a function of p, for k/u = 3 x 10~ 4 , 
calculated with the approximate propagator (dotted line) and with the obtained numerical 
propagator (solid line). 

the LO of our scheme, by using LPA' propagators. In the IR, we have seen that the 
result that we got for the anomalous dimension is 0.052 instead of a value close to 0.044 
that would have been obtained had we solved exactly eq. 1)16)1 . As a further test, we 
have recalculated I3 («; p) and J3 (p; «; p) using, instead of the LPA' propagators, the 
propagators © in which (p; p) is the 2-point function that has been obtained in this 
section by approximately solving eq. ()16)1. 

In fig. |U we plot the ratio of the function lf\n; p) (p ~ p/n, see. eq. f)A7)l ) calculated 
with the propagator obtained from the numerical integration of the flow equation divided 
by the function given by eq. 1)19)1. One can see that the smaller the value of k, the larger 
the difference, and that the main error is for values of p around the minimum p min of 
the effective potential (p m in goes from 1.8 to 3 as k runs from A to 0). Nevertheless, the 
difference stabilizes for small enough k and it never exceeds 4%. The right panel of fig. E] 
shows the comparison of the two curves in the worst situation, i.e., for small values of k, 

as a function of p: the curves are hardly distinguishable. 

(3) 

The same analysis is repeated for J 3 . Again, it is only for small values of k that 
the two functions differs. In the left panel of fig. El we display the ratio of the function 
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FIG. 5: Left: The ratio of the function Jg (p; k; p) calculated with the obtained numerical 
propagator and with the approximate LPA' propagators, as a function ofp/n, for k/u = 3 x 1CP 4 
and for different values of p. Right: The function J3 (p; k; p) / J3 (p = 0;k;p) calculated 
with the obtained numerical propagator compared to that calculated with the approximate 
propagators, as a function of p/k, for k/u = 3 x 10~ 4 and for different values of p. 

J3 calculated respectively with the obtained (numerator) and the approximate LPA' 
(denominator) propagators, for k/u = 3 x 10~ 4 , for various values of p. The difference can 
be large, but only in the region {p ^> k) where the function J3 itself is very small. In the 
region where the function is non negligeable, the difference between the two calculations 
never exceeds 5%. As was the case for the largest error occurs for values of p near the 
minumum of the potential. In the right panel of fig. we plot the two functions for the 
same values of k and p as in the left panel: the difference between the two calculations of 
is invisible on such a plot. 



V. CONCLUSIONS AND PERSPECTIVES 

We have demonstrated in this paper that the method proposed in Q| allows for 
concrete numerical applications. We have calculated the self-energy of the scalar model, 
at the LO of the approximation scheme, at criticality, at zero external field, in d = 3, 
and have obtained accurate results over the whole range of momenta. Already at this 
level of approximation the results obtained compare well with those of more elaborate 
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techniques. Worth emphasizing is the fact that the scaling behavior of the self-energy 
is accurately reproduced: not only do we get a reasonable estimate of the anomalous 
dimension, but the entire dependence of the self- energy on the momentum and the cut- 
off follows accurately the expected scaling behavior. To our knowledge, this is the first 
time that an approximate solution of the NPRG flow equations is constructed with these 
properties. 

In the present paper, whose main objective was to confirm the applicability of the 
method to a concrete calculation, we solved approximately the flow equation (|16|). How- 
ever, several tests suggest that this approximate solution differs in fact very little from 
the complete solution of Of course, a definite statement concerning the error made 
in the present calculation can only come from a comparison with the exact solution. This, 
we believe, is within reach. Similarly, work is in progress to test the convergence of the 
procedure by calculatin g th e next-to-leading order contribution. 



The method of ref. j.14] builds on our previous works on the same subject [15, 1 



The results presented in this paper indicate that it is both conceptually simpler, and 
numerically more accurate, than the method which we have developed in [Uj [ly]. It 
offers the possibility of applications to a variety of non-perturbative problems, where the 
knowledge of the momentum dependence of n-point functions is necessary. Even the 
approximate treatment presented in this paper could constitute an interesting starting 
point in situations where only a semi-quantitative description would be valuable. 
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APPENDIX A: THE p = SECTOR 

As discussed in the main text, our approximate solution of eq. (fTT)|) builds on the prior 
determination of quantities that are independent of momentum. These are calculated us- 
ing a variant of the derivative expansion that we describe in this appendix. The derivative 
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expansion is usually |7j formulated in terms of an ansatz for the running effective action 
r re [</>], including terms up to a given number of derivatives of the field. Its leading order, 
the so-called local potential approximation (LPA), assumes that the effective action has 
the form: 

T^LPArn f id ' ' 



[0] = * X l -d^ a d^ a + V K {p) | . (Al) 

where the derivative term is simply the one appearing in the classical action, and V K (p) 
is the effective potential. In the next-to-leading order (NLO), one assumes [lj 



^NLO\ 



d d X { ^Md^d^ + V K ( P ) } . (A2) 



An interesting improvement of the LPA, which we refer to as the LPA', is a simpli- 
fied version of the NLO that consists in ignoring the p-dependence of Z K (p), i.e., in 
chosing Z K = Z K (p ) where po is a given value of p, usually taken to be the run- 
ning minimum of the potential. In the LPA' one solves simultaneously the flow equa- 
tions for both the effective potential V K (p) (a partial differential equation in k and p) 
and for Z K . In this approximation, the inverse propagator takes the form of eq. (fTTj) : 
G^iq 2 ; 0) = Z K q 2 + V£(<f>) + R K (q), with V^{(j)) = d 2 V K /d<j) 2 . The LPA' allows for a non- 
trivial anomalous dimension, which is determined from the cut-off dependence of Z K (see 
eq. (j201) and ref. [l|). 

The procedure followed in this paper to determine the field renormalisation constant 
Z K differs slightly from that used in [l5(. This is because, as explained in sect. Ill, we 
need the calculation of Z K to be consistent with the approximate eq. ()16|) for the 2- 
point function. This is essential to get the proper scaling behavior of (p', p) at small 
momenta. We set: 

d^ K (p;p) 



Z K {p) = 1 



(A3) 

p=0 



dp 2 

where £ K (p; p) is defined in eq. (fi"5|) . Notice that, using eq. |T8|) . Z K (p ) = Z K . The flow 
equation obeyed by Z K (p) reads 

djf\p\ P ) 



nd K Z K (p) 



dp 2 



"V^^^.^,^ (A4) 



2 



p=0 

which follows immediately from eq. ()16j) for E K (p; p). Knowing the solution of this equation 
we can calculate r\ K from eqs. (}2T))) and (|T5|) . At this point, it is convenient to choose po = 0. 
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Then the expression of rj K that one deduces from eq. (jA4|) simplifies into: 

(A5) 



Since lf\p = 0) depends explicitly on Z K and r\ K (see eq. (|T9|) ). eq. (|A5|) is in fact a self- 
consistent equation for r\ K . The fact that a derivative of Z K (p) enters eq. (|A5|) demands 
the simultaneous resolution of eq. (IA4I) for small finite values of p. 

nn 

The solution of the LPA' is well documented in the literature (see e.g. |7, 20]). In 
practice, we work with dimensionless quantities. We set: 

v K (p) = K-/K- d V K (p), X (P) = ^A, (A6) 

with 

p = K d l Z K K 2 ~ d p, (A7) 

and Kd is given after eq. (|2Tj) . We solve the equation for the derivative of the potential 
with respect to p, i.e., w K (p) = dpV K (p), rather than that for the effective potential itself. 
This reads (from now on we stick to d = 3): 

K ^. = -(2-„.K + (1 -(l-W C^^ + 3< + 2M- 

V 5 / \{l + w K ) 2 {l + w K + 2pw'J\ 

(A8) 

where w' K = dpW K (p), w" = djw K (p). Eq. (jA8|) is solved starting from the initial condition 
at k = A: 

w K (p, k = A) = rh\ + g A p, (A9) 
where % and g\ are related to the parameters r and u of the classical action ((TJ) by 

™a = 9a = ^-jp (A10) 

and the parameter r is adjusted to be at criticality. Together with eq. (|A8|) . we solve the 
equation for x«(p > 0), which reads 

n~ / o~ (3< + 2M') 2 
«9«X« = + 1 + VkWk - 2 P / , , , - I 

(1 + < + p<) 4 

PXk1 5 J (l+< + 2p<)3 1 5 j (l + < + 2p<) 2 ' 1 J 
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FIG. 6: The dimensionless mass rh 2 K (f> = 0) = dpV K {p = 0) (left) and the anomalous dimension 
r\ K (right) as a function of k/u. These quantities were obtained by solving the LPA' equations, 
with u/A = 3.54 x and the parameter r adjusted to be at criticality. 

where x' K = dXn/dp. The initial condition is x«(p = 0) = 1 for all k, which follows from 
the definition of Z K , eq. (|TB*j) . Finally, for t] K we have simply: 

Vk = (i+w'J)ii X 'M/5 (A12) 

For the sake of illustration, we present in fig. [HI the LPA' solutions for rh 2 K (p = 0) 
(defined in eq. (J2B)) and r] K , as a function of k/u. The calculations have been done with 
u/A = 3.54 x 10 ~ 4 , but the curves are independent of this choice, provided u/A remains 
small. One can verify that the the crossover between the UV and IR regimes occurs 
around k c ~ u/10. The fixed point value of r\ K is rf = r] K ^ rs 0.05220. Fig. illustrates 
the p-dependence of the renormalization factor Z K (p) (see eq. f|A6j) ). This dependence is 
completely negligible when k ^ k c ~ w/10, and never exceeds 8%. 

APPENDIX B: THE FUNCTIONS lf ] (k; p) AND J 3 (3) (p;«;p) 

(2) (3) 

In this appendix we provide details about the functions J3 '(k; p) and J3 \p\K\p) 
calculated with the propagators ()17|) and ()22j) respectively. 

Consider first the function I^ 2 \k-; p), defined in eq. (fTTj) . and whose explicit expression 
is given in eq. The variation of 1$ (k; p) with k is dominated by the explicit linear 
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FIG. 7: The dimensionless function x{P) defined in eq. (|A6|) . for various values of k/u = 3 x 10 n , 
with n = 2 (circles), n = 1 (squares), n = (diamonds), n = — 1, • • • , —5 from bottom to top. 

k dependence and the ^-dependence of the renormalization factor Z K . The function 



(1 + ml(p)f V 5 

displayed in fig. |H1 illustrates the remaining dependence on k and p. 

Consider next J^ 3 \p; k; p), defined in eq. (jlTljl . Using the LPA' propagators of eqs. fPTj) 
and (|22|) one can calculate it analytically. One first makes the changes of variables p = 
p/n,q = qj k and cos 7 = p.q/p q, and then perform the integral over the remaining 
angular variables. One gets then 

1 1 

4 (p;«;p) = 72/o \2a \ ~2\2 / d ^ 2 / rf ( COS 7) 
KZZ(2n) 2 {l+ml) 2 J } _ x 

x (2 - 7/ + ryg 2 ) 

0(1 — q 2 — p 2 + 2qpcos r y) + (q 2 +p 2 — 2gj9cos7)©(g 2 +p 2 — 2gpcos7 — l)+m 2 

(B2) 

To perform the integral over cos 7 one needs to consider the various domains defined by 
the O functions. It is then convenient to separate the calculation in two different regions: 
p > 2 and p < 2, and this for the two possible signs of m 2 (see also (l5[). The calculation 
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FIG. 8: The function p)/n as a function of n/u, for different values of p. 



is then done by first integrating over cos 7; the remaining integration over q can be done 
by making first an integration by parts to get a rational function, that is then decomposed 
into simple fractions. One finally gets (the dependence on p is entirely contained in rh K (p) 
and is not written out explicitly): 

a) p > 2,m 2 K < 0. 

1 1 r_ 77/5 

2 V 3 

V V ( - 



2 + 



■f p 2 - 3m 2 K 



1 

+ 2p 



1 + \ + (p+ (1 - \ + \{p + V^K 2 



+ 



1 

2p 



7] 77 /_ 



log 



log 



■mt 



1 



A - 

KZl{2^{l+mlf {p* \3 



15 J t 



1 

15 



105 3 15 



p + 1 + \J—m\ 

p-l- \f ZI mi 
p+ 1 - v /Z ^f. 

+ ^(1/(P 6 )) 



(B3) 
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h)p<2,rhl <0. 



1 



+ 



^ 2 (27r) 2 (l+77* 2 ) 2 
1 



4 4 



25?? o 


+ 




1 


1 - 


V 


2p 




' 4 


1 


1 - 


V 


2^ 




' 4 



4 4?7 - . V-2 , 

1 + mi V 3 15 3 



8 



12 6) F 12(T 



log 
log 



3 i? 



P+ 1 + a/-?7?k 

1 + V /= ^I 



ml 



mi 



kZI{2*Y {I + ml)* {3(1 + ml) 

.2 , -2 



^ 2 + rj — 2m 2 . + rjih^ _ 3 2(1 + rj — 5m K + r/m K 



6(1 + m 2 J 3 
c)p>2,m2 >0. 

(3)' 1 1 



4 (p; 
i 



+- 



«^(2tt)2(1 + ^) 2 

^2 -2 



2 + 



15(1 + tb 2 ) 4 
5 
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p 4 + 0(p 5 ) 
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FIG. 9: The function Jg (k;p)/I^ (k) for p = (left) and p = 6 (right), as a function of (p/n), 
for different values of k/u. 

The function J3 (p;n;p)/J^ (p = 0; «; p) is displayed in fig. El for the two values 
p — 0, and p = 6. One sees that in both cases the p-dependence is concentrated in the 
region p ~ k: J3 (p; p) is independent of p when p ^ k, and it vanishes when p ^ a 
property that has also been exploited in and [la]. For p = 0, J^ 3 \k; p) / J^\k; p = 0) 
is essentially a function of p/n only. For p = 6 some residual dependence on k remains. 



APPENDIX C: ULTRAVIOLET BEHAVIOR OF THE SELF-ENERGY 

In this appendix we study the behavior of the self-energy £(p) for p 3> u. We show 
that the solution of eq. ([IB]) reproduces the result of 2-loop perturbation theory, namely 
£(p) = u 2 / (967r 2 ) log(p/w), albeit with a coefficient in front of the logarithmic that differs 
by 8%. 

Consider first the exact flow equation for the 2-point function, eq. (|5]L in vanishing 
external field (in this appendix p = throughout). At order 0-loop (indicated by the 
superscript [0]), this is simply: 

cUf)[°](p) = 0. (CI) 

This equation has the solution 

r^ 0] ( P )=p 2 , (C2) 

where we used the initial condition at k = A that one deduces from eq. (jT)), and adjusted 
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the bare mass r to be at criticality (S K=0 (p = 0; p = 0) = 0, yielding = 0). 

To go to 1-loop, one uses, in the r.h.s. of eq. J5J), the 0-loop expressions for both the 



propagator, G (k;p) = l/(p 2 + R K (p)), and the 4-point function T^°\p 1 ,p 2 ,P3.p4 : ) 
One gets 

dv^Uv) = -- f ddq ^(g) - u d f ddq 1 

KK lPJ ~ 2 J (2ttY (q* + R K (q)y 2 K J (2ir) d q 2 + R K (q) 

(C3) 

The integration is immediate; by imposing criticality and the initial condition at k — A, 
one obtains 

r<*, (p) . / + i /^{_L__^, (C4) 

giving a self-energy which is in fact independent of the momentum p. 

The 1-loop expression for the 4-point function, which will be needed shortly, is ob- 
tained similarly: 

d K T^\p, -p,l,-l) = 

u 2 J —^- d d K R R {q)Gl{K] q) {G {k; q) + G («; p + l + q) + G (k; p — l + q)}, (C5) 

which can be integrated easily to give 

T^(p,-p,l,-l) = 



u. 



a 2 



U — — / j^-^G (K-,q) {G (K-,q) + G (k-,p + I + q) + G (n-,p- I + q)} 



(C6) 



where we imposed the initial condition F^l A {p, —p,l, —I, p = 0) = u (the integrand in 
eq. (jC6|) should, for finite A, be subtracted from its value at k = A in order to satisfy this 
initial condition; the corresponding contribution, however, vanishes in the limit A — ► oo, 
and we assume here that A is large enough so that it can be neglected.) 

Going now to 2-loop, one puts in the r.h.s. of eq. © the 1-loop expressions of both 
the propagator and the 4-point functions. Since we are interested only in the momentum 
dependence of the 2-point function, we consider only the terms in the flow equation that 
depend on p, i.e., £ K (jo) = r^fjo) — T^^O) — p 2 . Since the momentum dependent terms 
originate entirely from the contribution of order u 2 in r( 4 )M, we can use Gq as propagator. 
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We have therefore 



«9 K £L 2l (p) = y J ^y d d K R K (l)G 2 ( K ;l) J -0-G o (K ; g)(G o («;p + / + g)-G o («;/ + g)). 

(C7) 

This expression can also be integrated to give 

sL 2] (p) = -y j W?S ^5 G? o( /c 5 Z )Go(«;?)(C?o(K;i, + Z + ?)-Go(«;/ + ff)). (C8) 

At this point, we need to deal with the fact that the 2- loop expression for the self- 
energy is IR divergent. And indeed when k — > at fixed p, the integral in eq. ()C8j) 
diverges. In order to go around this difficulty, we consider the derivative dpAT^ (p) 

8y} 2 Un\ 7/ 2 f d d l C d d n 

(C9) 

where p = p/|p| and R' K (q) = d q 2R K (q). The limit k — > can now be taken, and yields 

d^l (p) _u 2 f d d q 2q.p f d d l 1 1 

d\p\ 6 J (2n) d q A J (2vr) d I 2 (I + p - q) 2 ' 1 ' 

Performing the integral over I and those over cos# = p.q and |q|, one recovers the well 
known result ( in d — 3): 

d££Lo(p) 1 « 2 /' oc d|g| r jfl sin#cos# 

u 2 1 
96vr 2 |p|' 

Let us now turn to the perturbative limit of eq. (|16j). Note that, at both 0- and 1-loop 
orders, the predictions of eqs. (jHJ) and (|TKJ) for the self-energy coincide. A difference arises 
at 2-loop order since, at the LO of the approximation scheme, we should insert in eq. (|C7JI 
T^ l \p, —p, 0, 0) instead of ri 4 -"^(p, — p, I, —I) as we did in the exact calculation, where 
the expression of r^w is given in eq. (jC6|) . That is, the LO flow equation reads 

d^ ]LO ( P ) = Y j -^yd K Rn{l)Gl{^l) J -^- d G ( K ; q )(G ( K ;p + q)-G ( K ;q)). 

(C12) 



d6 

d\p\ 24 (27r)2 7 |g| J \p 2 + g 2 - 2\p\\q\ cosO) 1 / 2 

(Cll) 



26 



In contrast to what happens with eq. (jC7|) . here the integration over k can no longer be 
done analytically and we have to deal with a third integral. Let us call k' the variable of 
this integration, and integrate over t' = log(V/|p|). After making the changes of variables 
q — > \p\q and / — > \p\l, one obtains: 

4 2)[2]L °(p) = -yN 2(d - 3) 

f°° , f n n n^2 f A /" A 

x / dt 1 T> 

•/ log k/|p| 



^ 2/jr yd t rR K ,(l)G (K;l) J J^y G o(K; q)(G (K;p + q) - G (K;q)). 

(C13) 



Now, the derivative with respect to \p\ is very simple because, in d — 3, it only enters in 
the integration limit. One has: 

9sl 2)[2]LO (p) _ u 2 f d 3 l n mn2 , n /• A 



9jp| 2H7 (2vr) 3 



(C14) 



where t = log(«/|p|). In can be verified that the first term, i.e., that containing p, vanishes 
when k -C \p\. In this limit: 

u 2 1 

= AtA- (C15) 
9tt 4 |p| v ; 

Comparing eqs. (jCll|) and (|C15|) one sees that they both predict a logarithmic be- 
havior for the self-energy, the ratio of their respective coefficients being: 
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